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Unlike partially convective stars such as the Sun, fully convective stars do 
not possess a radiative core. Whether a star needs this core to generate a 


solar-like magnetic dynamo is still unclear. Recent studies suggest fully and 
partially convective stars exhibit very similar period—activity relationships, 
hinting that dynamos generated by stars with and without radiative cores 
hold similar properties. Here, using kinematic ages, we discover an abrupt 
change in the stellar spin-down law across the fully convective boundary. 
We found that fully convective stars exhibit a higher angular momentum 
loss rate, corresponding to a torque that is ~1.51 times higher for a given 
angular velocity than partially convective stars around the fully convective 
boundary. Because stellar-wind torques depend primarily on large-scale 
magnetic fields and mass-loss rates, both of which are suggested to be 
similar for partially and fully convective stars, the observed abrupt change 
in spin-down law suggests that the dynamos of partially and fully convective 
stars may be fundamentally different. 


The interiors of fully convective (FC) stars’ (M < -0.35 M,) and partially 
convective (PC) stars such as the Sun are fundamentally different, 
as fully convective stars do not possess a radiative core. How the 
stellar magnetic dynamo is affected by the existence of a radiative core 
and the role of the tachocline (the transition region between a star’s 
radiative core and its convective zone) is still unclear. On the obser- 
vation side, it is well known that the rotation periods and magnetic 
activities of stars are tightly correlated’, and recent observational 
studies have revealed that stars with different evolutionary stages? 
and masses* © exhibit similar period-activity relationships, which 
suggests that all stars obey a universal period-magnetic activity 
relationship and that a radiative core may not be a critical ingredi- 
ent for generating a Sun-like dynamo. On the theory side, various 
dynamo models suchas advection-dominated flux transport dynamo 
models’ ’ and turbulent dynamo models” can reproduce observations 
for both partially and fully convective stars to a certain extent, and 
thus this creates difficulties in understanding the types of dynamo 
operating in these stars. 


One way to break the degeneracy between theoretical predictions 
and observations is by understanding the time evolution of rotation 
periods of stars on either side of the fully convective boundary. Because 
adynamois ultimately responsible for generating the surface magnetic 
fields andthe magnetic activity that drives stellar winds, the amount of 
angular momentum carried away by stellar winds should be sensitive 
to the details of the dynamo processes. The evolution of rotation rate 
is sensitive to long-timescale trends in the average wind torque, which 
thus probes trends inthe dipole magnetic field strength and mass-loss 
rate, both of which should be tied to global dynamo relationships. 
Therefore, understanding the spin-down law across the fully convec- 
tive boundary could be the key to revealing the magnetic properties 
of stars and resolving this discrepancy. 

Theoretical rotation evolution models that are constrained by 
observed rotation period distributions have provided insight into 
the magnetic topology and angular momentumtransport instars" “. 
However, most of these works have focused on understanding FGK 
dwarfs, as both periods and ages for old M dwarfs are extremely difficult 
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Fig. 1| Spin-down of stars and stellar models near the fully convective 
boundary. Left, observational data showing the P,,,-age relationship for stars 

in the range 3,400 K < 7,4, < 3,500 K coloured by [Fe/H]. The marker styles show 
whether a star is determined to be FC (crosses) or PC (circles) using the stellar 
evolution model shown in the right-hand panel. Right, the fully convective 
boundary predicted by the STAREVOL”””’ stellar evolution model for stars with 
various Tand [Fe/H]. The shaded green area shows the temperature range of the 
stars in the left-hand panel, and the dotted red lines mark the typical metallicity 
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range in our sample as shown in the left-hand panel. The model suggests that 
stars with the same T, but lower [Fe/H] are fully convective whereas those with 
higher [Fe/H] are partially convective. Both PC and FC stars should exist between 
-3,400 K and 3,600 K. The bimodal distribution of rotation periods across stars 
of different metallicities seen in the left-hand panel suggests that the angular 
momentum loss rate changes abruptly across the fully convective boundary, in 
which the FC stars (top sequence; lower [Fe/H]) have a higher angular momentum 
loss rate compared with the PC stars (bottom sequence; higher [Fe/H]). 


to obtain. Old M dwarfs are faint and many rotate slowly (>25 days). 
This means that photometric data with high sensitivity and a long 
observational baseline are needed to measure their periods. Ages 
for old M dwarfs are also hard to infer as their observables change 
slowly with time, which creates challenges for age dating them with 
isochrone fitting. 

Various studies have provided hints on the spin-down of these 
low-mass M dwarfs towards older ages. Galactic kinematic and wide 
binaries ona relatively small sample of M dwarfs with periods obtained 
from MEarth’*° founda bimodality of fast and slowly rotating M dwarfs 
thatis difficult to explain with traditional models of angular momentum 
loss. However, the rotation periods measured for the 4 Gyr open cluster 
M67 (ref. 17) suggest that old M dwarfs do eventually converge onto 
a tight sequence. However, we still lack the sample size, especially at 
older ages, to constrain a spin-down law of these fully convective stars 
from observational data that can be used to test theoretical models. 

Arecent catalogue of rotation periods measured using the Zwicky 
Transient Facility (ZTF) has the sample size needed to understand 
the spin-down of fully convective stars. With the gyro-kinematic 
age-dating method”, we obtained kinematic ages for Kepler and ZTF 
stars with period measurements (see ‘Gyro-kinematic age sample’ for 
more details). We extended the age measurements for fully convective 
stars up to ~10 Gyr and detected an abrupt change in the spin-down law 
across the fully convective boundary. 


Results 

A double sequence near the fully convective boundary 

The left panel of Fig. 1 shows the rotation period (P,,,)-age relationship 
for stars with surface temperature (Tep) between 3,400 K and 3,500 K, 
coloured by metallicity taken from ref. 21. A bimodality emerges inthe 
P...-age relationship in this temperature range with metal-poor stars 
spinning down quicker than metal-rich stars. 

Stars with similar temperatures can have either fully or partially 
convective interiors depending on their metallicity. Figure 1 (right-hand 
plot) shows the fully convective boundary (red line) for stars with dif- 
ferent [Fe/H] and temperatures predicted from the STAREVOL stellar 
evolution model” (for details, see ‘Stellar evolution model’). We 
define the fully convective boundary to be the mass at which stars are 
able to maintain a radiative core after 10 Gyr. Below this threshold, 


there is a small range of mass where stars undergo the so-called ‘kiss- 
ing instability’ that usually brings them back toa fully convective state 
after a few cycles”. The model predicts the range of masses for stars 
that go through ‘kissing instability’ to be -0.03 M,, and hence the 
width of the fully convective boundary to be ~0.03 M,. Note that we 
selected the surface temperature of each model at about 5 Gyr. The 
model suggests that, near the fully convective boundary, stars with the 
same temperature but lower [Fe/H] are FC whereas those with higher 
[Fe/H] remain PC, and that both fully and partially convective stars 
exist between ~3,400 K and 3,600 K. This means that, given a sample 
of stars with a wide range of [Fe/H] (since the ZTF survey covers the 
entire northern sky, this statement should be satisfied), we can com- 
pare the P,,,-age relationship for both fully and partially convective 
stars with similar masses simultaneously by selecting stars ina narrow 
range of temperatures near the fully convective boundary. Because the 
upper sequence in the left-hand plot of Fig. 1 is more metal poor than 
the lower sequence, we interpret this to mean that the upper sequence 
consists of FC stars and the lower sequence consists of PC stars. It is 
also worth pointing out that metallicity measurements are correlated 
to age, mass and rotation. This means that the separation of the two 
sequences could be due to factors other than [Fe/H]. However, as shown 
in the next section, the observable feature in the colour-magnitude 
diagram (CMD) that is associated with the fully convective boundary 
and is independent of the metallicity measurements’ clearly sepa- 
rated the two sequences. 


The bimodality of spin-down laws 

To further understand whether this bimodality is actually caused by an 
abrupt change in the spin-down law of stars across the fully convective 
boundary, we separated the fully and partially convective stars based 
ontheir absolute Gaia magnitude, M, and Gaia BP-RP colour measure- 
ments using Jao’s gap”. We used the Jao gap to separate the stars because 
the measurements for M, and Gaia BP (G,,)—Gaia RP (Gpp) colour are 
reliable. Jao’s gap is an under-density inthe CMD near the fully convec- 
tive boundary discovered using stars within 200 pc of the Sun from 
Gaia DR2 (ref. 26). This gap can be approximated by a line connecting 
(Mg, Ger — Grp) = (10.09 mag, 2.35 mag) and (Mg, Gsp — Grp) (10.24 mag, 
2.55 mag) and is thought to be caused by structural instabilities due 
to non-equilibrium fusion of °He (refs. 24,27). In the rest of the paper, 


Nature Astronomy 


Article 


https://doi.org/10.1038/s41550-023-02126-2 


3,200 K < Ty < 3,300 K 


3,300 K < T, < 3,400 K 


P ot (days) 


Partially 
convective 


84.83 days 


65.52 days 


7517 days 


Po (days) 


51.03 days 


7517 days 


L 46.21 days 


Gyro-kinematic age (Gyr) 


Fig. 2 | Spin-down sequences for fully and partially convective stars in narrow 
temperature bins. P,,.-age relationships for stars in four different narrow Ter 
bins, where the fully convective stars (red) and the partially convective stars 
(black) are separated using Jao’s gap”. The five black lines are Skumanich spin- 
down laws” (P,,.« age°>) to guide the eyes. These lines show that, at a given age, 
the fully convective stars have a spin period that is ~1.5 times larger than the 
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partially convective stars. The marker styles indicate whether the star is fully 
convective (empty circle) or partially convective (filled circle) based on the 
STAREVOL stellar evolution model. The normalized histograms (bottom axis 
ranges from 0 to 0.05) onthe right of each subplot are the period distributions 
for fully convective stars (red) and partially convective stars (black) that have 
gyro-kinematic ages >2 Gyr. 


we use this line to roughly separate the fully and partially convective 
stars; stars lying above this line in the CMD are more likely to be partially 
convective and those below are more likely to be fully convective. 

We plotted the P,,,-age relationship for FC (red) and PC (black) 
stars in narrow temperature bins of 100 K, between 3,200 Kand 3,600 K 
(Fig. 2). The normalized histograms on the right show the bimodal 
period distributions of these stars with gyro-kinematic age >2 Gyr and 
the lines mark the bins with the highest normalized number density. It 
is notable that the double sequence mostly exists for stars in the range 
3,300 K< Tee < 3,500 K: this is slightly lower than, yet very close to, the 
Terange that the stellar evolution model predicts to contain both fully 
and partially convective stars if a sample contains stars with a wide 
range of [Fe/H] (Fig. 1, right-hand plot). More interestingly, stars below 
3,300 K (fully convective) follow the top sequence and those above 
3,500 K (partially convective) follow the bottom sequence, which 
further suggests that the spin-down laws of fully and partially convec- 
tive stars are bimodal. However, it is worth pointing out that this gap 
does not provide a clean division between partially and fully convective 
stars as they can oscillate between this gap while transitioning between 
being partially and fully convective”. 

The abrupt change in the rotational evolution between fully and 
partially convective stars means that the angular momentum loss rates 
also exhibit an abrupt change between stars with and without a radia- 
tive core. As shown by the detailed calculation in ‘Torques for FC and 
PC stars spinning at the same rate’, a fully convective star experiences a 
higher spin-down torque by a factor of ~1.51 than a partially convective 
star with the same rotation period. Stellar-wind theory then predicts 


that, at a given rotation period, fully convective stars should have 
dipole fields that are ~1.26 times stronger or mass-loss rates that are 
~1.44 times higher (or some combination of both). 


The mass dependence of the spin-down laws 

With gyro-kinematic ages, we can directly examine the rotation period 
distributions and spin-evolution isochrones for partially and fully con- 
vective stars. We fitted a Markov chain Monte Carlo model using 
emcee” by maximizing the log likelihood, calculated by 
—0.5 >) (log oP rot — log ,oPrit) where Pais given by 


b 
M, 1 days 
= cy 70.5 * < 
Par = L097 (sax) Gyr for M, < Mfc (1) 
b 
M, 2 days 
= 0.5 * 
Pre = 1021 (sasiz) Gyr for M, > Mec (2) 


where M;cis the mass of the fully convective boundary (~0.35 M, based 
on stellar evolution models), Tis the age of the star in units of billion 
years and M. is the mass of the star in units of M,. The best-fitted param- 


eters are b, = —0.770°77 , , by = —0.9380448  , Cı = 1.4489.958 and 
Cy = 1.395000 Mpc = 0.331900 Mo. The corner plot is shown in Sup- 


plementary Fig. 1. Itis worth pointing out that the uncertainty for the 
mass for the fully convective boundary is similar to the range of masses 
that gothrough kissing instability predicted from STAREVOL (-0.03 M,). 
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Fig. 3 | Period-mass diagrams plotted with Markov chain Monte Carlo best- 
fit models. log,.(P,..)-mass distribution for the entire sample (left), partially 
convective stars (middle) and fully convective stars (right) coloured by gyro- 
kinematic ages. The green solid lines are isochrones every 1 Gyr apart from 1 to 


10 Gyr, derived from equations (1) and (2). The best-fit model suggests that the 
spin-down law for fully convective stars is slightly more mass dependent than for 
partially convective stars. 


Figure 3 shows the log,o(P,..)-mass diagram for all the stars (left), 
partially convective stars (middle) and fully convective stars (right) 
coloured by gyro-kinematic ages. The masses are calculated using 
STAREVOL with metallicity”. The green lines show the predicted isochro- 
nes every 1 Gyr apart from 0 to 10 Gyr using the model described above. 

The intermediate period gap, an observed dearth of stellar rota- 
tion periods in the temperature-period diagram at ~20 days for G 
dwarfs and up to ~30 days for early-M dwarfs, is believed to be caused 
by stalling of stars due to core-envelope coupling”. Separating the 
stars into fully and partially convective shows that this gap appears only 
in the partially convective stars, which further supports the hypoth- 
esis that this period gap is formed from the redistribution of angular 
momentum between astar’s radiative core and its convective envelope, 
as fully convective stars do not possess a radiative core. 

Based on the model fit, the fully convective boundary is at 
-0.331 M,. The powers of the mass-dependent part of the spin- 
down law agree within 1o for fully convective stars and partially convec- 
tive stars. At the fully convective boundary, fully convective 
stars spin-down faster than partially convective stars by 
[10% (Mrc/0.35M,)"|/[10°(Mrc/0.35My)"2] = 1.23. This gives a 1.51 
times higher angular momentum loss rate for fully convective stars, 
which corresponds to a magnetic field strength that is 1.26 higher, a 
mass-loss rate that is 1.44 higher or a combination of both”. For more 
details, see ‘Torques for FC and PC stars spinning at the same rate’. 


Discussion 
The bimodality of spin-down laws for fully and partially convective 
stars means that they have fundamentally different spin-down laws 
and thus different angular momentum loss rates. Observational data 
suggests that fully convective stars lose angular momentum -1.48 times 
faster than partially convective stars at a given angular velocity (see 
‘Torques for FC and PC stars spinning at the same rate’ for more details 
about this calculation). Because stellar spin-down and winds are direct 
consequences of the stellar dynamo, this suggests that the dynamos of 
fully convective and partially convective stars are also fundamentally 
different. This result is pronounced observational evidence that the 
dynamos of fully convective and partially convective stars operate 
differently. However, the exact operational difference between their 
dynamos is not clear, as dynamos and the mechanisms that convert 
magnetic energy into the heating that drives stellar winds and angular 
momentum loss are still not well understood. 

Typically, stellar-wind theory indicates that the angular momen- 
tum loss rate should be directly correlated with the mass-loss rate, and 


the strength and geometry of the magnetic fields*’*. Observations”? 
suggest that the mass-loss rate for fully convective M dwarfs is similar 
or smaller than that of the Sun based on the ultraviolet spectra of stellar 
H 1Lya lines from the Hubble Space Telescope. Activity indicators such 
as X-ray” and Ha (for example)? also suggest that fully and partially 
convective stars exhibit similar Ro-activity relationships, which indi- 
cates similar magnetic field strengths. However, our result suggests an 
abrupt change in the wind torque of only ~50% (shown in ‘Torques for 
FC and PC stars spinning at the same rate’) and it is not yet clear what 
may be different in the stellar winds to produce that change and thus 
whether or not sucha change should be visible above the scatter in the 
observed Rossby number (Ro)-activity relationships. 

If the magnetic field strengths are indeed similar, by combining 
previous studies and this work, we speculate that the differences inthe 
dynamos of fully convective and partially convective stars exist in their 
magnetic morphology, which causes fully convective stars to generate 
magnetic dipoles that are stronger but of similar higher-order magnetic 
fields compared with partially convective stars. Angular momentum 
loss through magnetic winds (stellar spin-down) is likely to be driven by 
the escape of open field lines produced by the magnetic dipole**. The 
abrupt change in spin-down law across the fully convective boundary 
discovered in this paper suggests that fully convective stars are likely 
to have a more poloidal magnetic field configuration than that of 
partially convective stars. 

That fully convective stars might exhibit different dynamos than 
partially convective stars is also hinted at when looking at fast-rotating 
M dwarfs”. These fast-rotating, fully convective M dwarfs have either 
predominantly dipolar, axisymmetric global fields or multipolar, 
non-axisymmetric dipolar fields. 


Methods 

Gyro-kinematic age sample 

We determined gyro-kinematic ages following the procedure described 
in ref. 20, where the vertical velocity dispersion for each star was cal- 
culated from vertical velocities of stars that are similar in temperature 
(Tas calculated from G,p — Grp measurements using a polynomial fit 
taken from ref. 30; Gsp — Grp de-reddened using dustmap*”), rotation 
periods (P,,,), absolute Gaia magnitude (M,; extinction-corrected using 
dustmap) and Rossby number (Ro) to the targeted star. The turnover 
times used to calculate the Rossby number were taken from See et al. 
(manuscript in preparation) and have the form 


Tez = 0.388 x 10% ent@2Ter*+% for Taye < 3,480K (3) 
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Tez = 0.388 x 10% TertbzTartbsTertba for Tope > 3, 480K. (4) 
This form was found by fitting polynomials to log Tez and Teg from the 
stellar structure models” at zero age main sequence. The values of the 
coefficients are a, = 6.5211 x 107, a, = -4.0036 x 10”, a, = 8.6823, 
b, =-2.5190 x 10, b, = 3.7361 x 10°, b, = -1.8557 x 10 and b, = 32.5951. 

Wethen converted the velocity dispersions into stellar ages using an 
age-velocity-dispersion relationship”. Post-main-sequence stars were 
cut out by only selecting stars with M, > 4.2 mag and equal-mass binaries 
were excluded by fitting a sixth-order polynomial to the main-sequence 
stars in M,-Tog space, which moved the fitted polynomial along the M, 
axis so that it lay directly below the equal-mass binary sequence and 
removed stars with M, greater than the modified polynomial. 

The dataset used in this work is from Lu et al.” and combined the 
-30,000 stars in ref. 20 from Kepler*’ and the -55,000 stars with period 
measurements from ref. 18 and Lu et al.” from ZTF”. The vertical veloci- 
ties for the ZTF stars were obtained using radial velocity measurements 
from Gaia DR3 (ref. 42). We did this by transforming from the Solar 
System barycentric International Celestial Reference System refer- 
ence frame to Galactocentric Cartesian and cylindrical coordinates 
using astropy****. The bin size was (Tor, 10g10(P;ot), RO, Mc) = (177.8 K, 
0.15, 0.15, 0.2 mag), optimized following the methodology outlined 
in ref. 20 using clusters of stars ranging from 0.6 to 4 Gyr (refs. 17,30). 


Stellar evolution model 

The stellar evolution model for this work was computed using 
STAREVOL””?. We used a refined grid of standard models for masses 
between 0.3 and 0.4 solar masses by steps of 0.01 for eight metallici- 
ties between [Fe/H] = -1 and +0.3. Abundances were taken from ref. 45. 
We used an analytical surface atmospheric fit*’ and a solar-calibrated 
mixing length parameter a = 2.11. 


Eliminating possible systematic causes 
To eliminate the possibility of systematic effects or biases artificially 
producing this result, we visually examined 100 random ZTF light 
curves for the fully convective stars between 3,300 K and 3,500 K 
(where the double spin-down sequence exists) with measured periods 
>50 days. We found no systematic effect inthe measured periods that 
could cause all the fully convective stars to be on a period harmonic. 
Moreover, rotation periods in the spin-down sequence for the fully 
convective stars are not period harmonics (integer multiples) of those 
for the partially convective stars (Fig. 1). Itis also unlikely that a bias or 
systematic effect in rotation period measurements would conspire to 
create a double sequence only around the fully convective boundary 
(Fig. 2) predicted by the stellar evolution model (Fig. 1) and that the 
CMD gap is able to nicely separate the two sequences (Figs. 2 and 3). 
To eliminate the possibility of a systematic effect that could be 
caused by combining the two data sets from Kepler and ZTF, we per- 
formed the same tests with just the ZTF sample with radial velocity 
measurements from Gaia DR3 and found no notable changes in the 
results presented in this paper. This is expected, as Kepler mainly tar- 
geted solar-like stars and thus provided only 198 stars in the tempera- 
ture range we are interested in (3,200-3,600 K) compared with 7,809 
stars from ZTF. As the number of ZTF stars is more than a magnitude 
higher than that of Kepler, it is expected that excluding Kepler stars 
would not affect our results in this study. 


Limitations of the gyro-kinematic age-dating method 

The gyro-kinematic age-dating technique assumes that stars with 
similar parameters (effective temperature, rotation period, Rossby 
number and absolute magnitude) are approximately the same age. In 
general, it is expected that this assumption should hold. However, if 
partially convective and fully convective stars have different braking 
laws, this assumption may be broken for stars near the fully convec- 
tive boundary. Stars on either side of the boundary, with otherwise 


similar properties, could have quite different ages. It is therefore worth 
examining the behaviour of the gyro-kinematic age-dating method at 
the fully convective boundary in more detail, given that our results rest 
upon this technique. 

We considered a star near the fully convective boundary; whether 
it is on the fully or partially convective side of the boundary does not 
matter. To estimate a gyro-kinematic age for this star, we calculated the 
vertical velocity dispersion of all the stars contained ina bin in param- 
eter space centred on the star in question and then converted this to an 
age estimate using an age-velocity-dispersion relationship. Because 
we are near the fully convective boundary, this bin could contain both 
fully convective and partially convective stars. If fully convective and 
partially convective stars do indeed obey different braking laws, then 
the fully convective and partially convective stars are likely to have 
different ages even though they are contained within the same bin. 
The estimated age for the star will therefore be an average age of the 
fully convective and partially convective populations. This would have 
the effect of smearing out the sequences in Figs. 1 (left) and 2. The fact 
that we still see two sequences in these figures, even after this smear- 
ing out, suggests that our results are robust to this limitation of the 
gyro-kinematic age-dating technique. Indeed, in reality, the sequences 
could be even more well defined than shown in this work. 


Torques for FC and PC stars spinning at the same rate 

Inthis section, we determine how much larger the braking torque acting 
on FC stars is compared with PC stars by considering their observed 
rotation evolution. At late ages, it is thought that stellar rotation peri- 
ods increase with the square root of age, which is the well-known Sku- 
manich relationship, 


P(t) = at, (5) 


or, interms of angular frequency, 


Ab = Myos, (6) 


Here, ais a constant of proportionality that we will empirically deter- 
mine later. The rotation evolution of low-mass stars is governed by the 
angular momentum equation 

(7) 


dQ _T 
7 


do _T_ dt 
dt I 


dt’ 


where Tis the spin-down torque and /is the moment of inertia. On the 
main sequence, the d//dt term is approximately zero because a star’s 
stellar structure does not appreciable change during this phase of evolu- 
tion. Therefore, by differentiating equation (6), substituting into equa- 
tion (7) and rearranging, one finds that the spin-down torque is given by 


(8) 


The ratio of the torques acting on FC and PC stars, assuming the same 
angular frequency, is therefore 

Tre = (aye (9) 
Tec = \@pc/ Ipc 


For stars at the fully convective boundary (M. = 0.331 Mo), Qc/Qpc = 1.23 
(see main text, derived from equations (1) and (2)). As a result, if we 
consider the case where /,,//p- = 1, we find that the ratio of torques is 
T;c/Tpc 7 1.23” = 1.51. For this to be the case, stellar-wind theory" predicts 
that fully convective stars should have dipole field strengths that are 
larger bya factor of ~1.26, mass-loss rates that are ~1.44 times larger or 
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some combination of both of these factors. Further observations are 
needed to determine which of these scenarios is true. 


Data availability 


The data used are available in Supplementary Data 1. 
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